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Abstract 

Particle-in-cell simulations were performed by using 
PARMELA to characterise an electron injector with a 
booster linac for the AWAKE project in order to provide 
the baseline specifications required by the plasma wakefield 
experiments. Tolerances and errors were investigated. A 
3 GHz travelling wave structure designed by using CST code. 
Particles were tracked by using the field maps acquired from 
these electromagnetic simulations. These results are pre¬ 
sented in comparison with the generic accelerating structure 
model within PARMELA. 

INTRODUCTION 

The AWAKE project is a proton driven plasma wakefield 
acceleration experiment by utilising the driver proton beam 
from CERN’s SPS injector and a custom photo injector for 
the witness (trailing) electron beam m- 

Eigure shows the layout of the injector consisting of 
an S-band standing wave RE gun (SW), previously used in 
CERN’s PHIN photo injector Q, followed by laser optics to 
direct the laser onto the cathode. Beamline continues with 
beam current and position monitors and a pepper pot emit- 
tance measurement system suitable for a beam subject to 
space charge effects. A new S-band booster section (acceler¬ 
ating travelling wave structure, ATS) was designed by using 
CST suite 0 and introduced in the model to boost the beam 
energy to a region adjustable between 16 — 20 MeV. ATS is 
followed by a quadrupole triplet and a downstream screen 
to perform quadrupole scans for emittance measurement. 

The tracking studies have been performed by using 
PARMELA Q. For ATS the model provided by PARMELA 
and field maps extracted from CST were used and compared. 
Implementation of CST maps into PARMELA are discussed 
in the following sections. 

BOOSTER LINAC 

An S-band booster linac, ATS, was designed as a trav¬ 
elling wave structure with constant gradient of 15 MV/m 
through the entire structure (Figj^ . It consists of 30 cells 
with 120° phase advance and varying radii matched to 1 /im 
precision. ATS was optimised for low reflection coefficient 
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of about 2.5%. The multipole terms |[^ due to transverse 
RF-kicks are 9.4241 x 10“^ mT, 7.8418 x 10“^ mT/m and 
4.9 X 10“^mT/m^, respectively, from dipole to sextuple 
terms. 



Figure 2: A snapshot from CST design of the S-band trav¬ 
elling wave structure as the energy booster linac for the 
electron injector. 

USING CST FIELD MAPS IN PARMELA 

In order to use CST field maps in PARMELA tracking 
simulations, a MATLAB script was prepared to for¬ 
mat the standard CST field maps into the form required 
by PARMELA . The information that PARMELA requires 
on each line of a field map can be found in Table IV-2 from 
the program manual Q . 

In PARMELA, two field maps must be provided for a trav¬ 
elling wave structure; one produced with Neumann boundary 
condition (cosine map) and the other with Dirichlet boundary 
condition (sine map). These fields which are shifted in phase 
by 90° are fed into PARMELA by using the TRWCFIELD 
command. A single TRWAVE line is used to represent the 
entire ATS including the bore tubes with lengths equal to 
a cell length at each end of ATS to account for the fringe 
fields. 

BASELINE DESIGN 

In order to maintain the balance between the emittance 
growth and the energy spread within ATS, a 10° off-crest 
phase was chosen for the baseline design by using the CST 
maps of the current design while the SW is optimised for 
the highest energy. In addition, the same beam dynamics 
conditions can be met with on-crest RE phase in the case of 
the PARMELA model for a travelling wave structure. 

The space charge component of the emittance is com¬ 
pensated using the field produced by two solenoids located 
around the SW. Slightly different working points to satisfy 
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Figure 1: The layout of the photo injector and the booster linac section to provide a witness beam for the AWAKE plasma 
Wakefield acceleration experiment. 
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Figure 3: a) Transverse normalised emittance and b) the 
beam size as a function of the distance from the cathode. 
(Solid and dashed lines are x and y axes, respectively.) 

the constant envelope condition were determined for the field 
distributions from PARMELA and CST models. 

The beam dynamics specifications are presented in Table 
[^in comparison with these two. 

RF PHASE ERRORS 



Figure 4: Evolution of the beam energy spread as a function 
of the distance from the cathode. 


Table 1: Required specification and values produced with 
different models. 


Parameter 

Required 

PARMELA Map 

CST Map 

E (MeV) 

16 

16 

16 

AE/E (%) 

0.5 

1.3 

1.1 

ex,y (mm mrad) 

2 

1.7,1.8 

2.3, 2.5 

Px,y (m) 

5 

4,4 

0.8, 0.8 


0 

-0.2,-0.2 

-1.4,-1.2 


the RF jitter as the error source. Resulting variations were 
presented in Table Especially time of flight errors are 
crucial to determine as they affect the synchronisation in the 
entrance of the plasma cell. 


Table 2: Implications of the phase errors on various beam 
dynamics observable. 


Parameter 

1(7 Error 

Error/Degree of the phase jitter 

E (keV) 

1.4 

5 

AE/E (ppm) 

28 

93 

ex,n (nm-rad) 

2 

7 

Time of flight (fs) 

57 

190 


Variations in the beam dynamics due to the RF phase er¬ 
rors induced on SW and ATS were investigated. As they are 
planned to be fed by the same klystron, the same phase error 
can be assigned to both structures. A typical phase error 
of 300 fs (Icr) was induced over a 500 samples assuming 


EMITTANCE DIAGNOSTICS UNDER 
SPACE CHARGE 

For a photo injector, the dominant limitation for the emit¬ 
tance diagnostics is the effect of space charge. Figure 
















































shows the ratio of two defocusing terms from the envelope 
equation; the space charge and the outward pressure due 
to beam emittance. Different curves correspond to differ¬ 
ent bunch charges and lengths which are considered within 
the operational range of the injector. Although some cases 
are less space charge dominated (> 1) than others, in all 
cases considered, space charge is the dominating component 
of the total beam emittance. In the pre-booster region, a 



Figure 5: The ratio of space charge and emittance term of 
the beam envelope as a function of the distance from the 
cathode. 
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7: Beam size, divergence and the correlation terms 
range of working points of the injector. 


developed to use sparse data from a quad scan to measure the 
rms emittance of space charge dominated beams Thus, 
the lower limit of measurable divergence values is 0.3 mrad 
for the injector. Figure [ 7 ] shows the beam size, divergence 
and the correlations terms for a range of working points for 
the injector. As seen in the figure, the divergence values 
range from 0.5 and upwards implying that the ODR-DFRI 
technique could have been implemented for this system as a 
complementary technique to the conventional quadrupole 
scan in case of more severe space charge limitations. 


pepper-pot measurement system will be implemented as a 
standard technique to measure the emittance of space charge 
dominated beams. It is currently being designed at the Cock¬ 
croft Institute. For the post-booster region quadrupole scan 
technique is planned. The feasibility and reliability of this 
technique was assessed by simulations. Figure [^compares 
the emittance values calculated from a simulated quadrupole 
scan and the rms emittance value at the same location. The 
results differ 2% implying that emittance can be reliably 
measured within measurement error by using quadrupole 
scan in the post-booster region. 



Figure 6: A simulated quadrupole scan result in the post¬ 
booster region. 

We have also evaluated the Optical Diffraction Radia¬ 
tion Dielectric Foil Radiation Interference (ODR-DFRI) to 
measure the rms emittance. This technique which is a modi¬ 
fied version of the Optical Transition Radiation Interference 
(OTRI) technique extends the latter to lower beam en¬ 
ergies and has the advantage of single shot operation. The 
technique can be used in conjunction with a new algorithm 


CONCLUSIONS AND OUTLOOK 

The tracking simulations for the injector were performed 
by using the field models from both PARMELA and CST. 
Slight discrepancies in the results are under investigation. 
Errors due to phase jitter were studied. Emittance diagnos¬ 
tics under space charge effect is studied for both low and 
higher energies. 
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